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NOMENCLATURE 


u 

t 

a 

d 

M(t) 

B(t) 

I 


velocity  potential  of  flow 
acceleration  potential 
free  stream  velocity 
time 

normalized  cavity  pressure 

fluid  density 

slope  of  foil 

y  Intercept  of  foil 

semichord  of  foil 

position  of  rear  end  of  cavity 


a 


A 


B 

f 

m 


(0 

f_,m_,a>_ 
o*  o’  o 

7 

0 


force  on  foil 
moment  of  foil 

cavity  length/chord  length 

steady  parts  of  f,m,a) 
frequency  of  oscillation  of  foil 

reduced  frequency  <^) 

Bessel  functions  (see  [3]  page  172  ) 


Iv 
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? 

r 

Theodorsen's  function 

5k 

6 

net  normalized  unsteady  pressure  on  foil 

f* 

normalized  unsteady  force 

* 

m 

normalized  unsteady  moment 

* 

a> 

normalized  imsteady  cavity  length 

f*  f* 
s’  c 

real  and  Imaginary  parts  of  f 

*  * 

S’  c 

real  and  Imaginary  parts  of  m 

*  * 

“s’  “c 

real  and  Imaginary  parts  of  co 

^0 

steady  part  of  M(t) 

Ml,  Bi 

anq)lltudes  of  unsteady  parts  of  M(t),  B(t) 

The  following  quantities  have  complex  definitions 

given  In  (1] 

pages  4  and  7 

^k’4 

Appendix  VI 

Qi’Q2 

page  30 

page  34 

A*  A* 

‘Hc,c’  \b 

page  32  (see  also  page  22) 

“c’  “s 

page  22 

•^k,c’  ‘^k,s 

page  17 

^k,c’  ^k,s 

page  22 

^a* 

C  Ja,s-  ^a.o- 

V 
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1. 


1.  SUMMARY 

This  report  contains  the  results  of  calculations  of 
the  unsteady  force,  moment,  and  cavity  length  of  a  partially 
cavltated  hydrofoil  subject  to  unit  heave  or  pitch  oscillations. 
These  can  be  combined  with  arbitrary  amplitude  and  phase. 

Included  also  are  analytic  descriptions  of  the  zero-cavity- length 
limit  and  of  the  zero -frequency  limit.  Numerical  results  for 
these  limiting  cases  are  also  Included. 
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2.  INTRODUCTION 

A  previous  repoirt^^^  presents  the  theory  for  the  un° 
steady  motion  of  a  partially  cavltated  hydrofoil.  In  order  to 
carry  out  a  hydroelastic  analysis  and  also  to  Indicate  directions 
for  further  experimental  and  theoretical  work,  numerical  calcula¬ 
tions  had  to  be  made.  The  results  of  these  calculations  form 
the  main  body  of  this  report. 

For  certain  limiting  cases,  namely  the  cavity  length 
approaching  zero  and  the  reduced  frequency  approaching  zero, 
we  obtain  simple  expressions  for  the  solutions.  For  the  cavity 
length  going  to  zero,  we  show  that  the  limiting  result  agrees 
with  Theodorsen's  standard  result 

The  limiting  result  for  zero  reduced  frequency  was 
obtained  by  a  quasi-steady  analysis  of  the  steady  state  solution. 
The  results  agree  graphically  with  the  limit  of  the  results  of 
the  above-mentioned  calculations. 

A  short  summary  of  [1]  Is  Included.  This  makes  It 
possible  to  use  the  results  presented  herein  without  reading 
all  of  [1]. 
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3. 


3.  SIMKLRY  or  T»C-153-SR»1 

In  [1]  the  problem  under  study  was  the  determination 
of  the  flow  around  a  partially  cavitated  hydrofoil.  The  foil 
was  asstimed  to  be  a  flat  plate  of  Infinite  span  and  immersed  in 
an  infinite  fluid.  The  free  stream  velocity  and  the  cavitation 
number  were  assumed  constant.  The  foil  was  assumed  to  be  subject 
to  unsteady  motion,  specifically,  simple  harmonic  heaving  and/or 
pitching,  and  the  particular  question  raised  was  the  determina¬ 
tion  of  the  unsteady  terms  for  the  flow  in  these  cases. 


Tioman's  general  solution  in  terms  of  the  accelera¬ 
tion  potential  defined  by  ^  +  Ui^  where  i  is  the  velocity 

potential,  is 
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^  -  a  +  2 

k-1  K 


where  ^  are  explicitly  given*  (Appendix  VI),  and  are 
defined  by  the  motion  of  the  foil  (page  7),  while  and  C2 
are  (along  with:  tile  length  of  the  cavity)  to  be  determined. 
The  four  possible  boundary  conditions  to  determine  the  three 
unknown  functions  are  stated,  (page  6) .  The  vertical  com¬ 
ponent  of  the  velocity  1^,  Is  derived  In  terms  of  the 
conjugate  functions  of  to  be: 


C 


The  steady  state  solution  Is  then  derived  (pages  9-12) 
with  the  result 


Cj  -  -  2^  (3a*  +  1  +  SaVa^+l) 


with  a  (defined  by  a  * 


)  given  by: 


(^/?+I  +  a)^  . 


A  second  linearisation  Is  carried  out  and  the  problem 
Is  reduced  to  solving  six  linear  equations  In  six  unknowns 
(pages  13-24),  where  the  six  equations  are  given  on  page  23. 

The  unsteady  lift  and  moment  are  derived  In  terms  of  the  unknown 

*  All^^pendlx  and  page  number  references  In  this  section  are 
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constants  (pages  25-28)  and  are  given  on  page  28,  along  with 
the  steady  lift  and  moment. 

In  pages  29-35  we  rewrite  the  equations  In  a  non- 
dlmenslonal  form  suitable  for  machine  computation.  We  define 
solutions  for  unit  heave  and  unit  pitch  as  functions  of 
(Djj  ■  cavity  length/chord  length  and  SI  -  reduced  frequency,  and 
show  how  to  combine  these  solutions  for  arbitrary  pitch  and 
heave. 
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4.  ZERO  CAVITY-LENGTH  LIMIT 


For  analytical  purposes  It  Is  convenient  to  represent 
the  set  of  six  real  linear  equations  by  three  complex  linear 
equations.  Let  Q*  ■  Q2  +  IQj^,  v*  ■  +  Iv^,  and  let 

ic 

q  -  qy  +  Iqj^  where 


where 

<  -  +  <c 

«k  -  •'k.s  +  ^ 

-  -^(•’k.c  +  ^  +  l»k,s  +  ^  -  i^k<")» 


+  «k,.>  -  -  K 
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7. 


Therefore 


i  ZdT 
3 


In  this  notation 
Q  V  -  Iq 


Let  rj^  - 


To  determine  the  asymptotic  properties  of  v  ,  It  Is 

necessary  first  to  determine  the  asymptotic  properties  of  Q  , 

r^,  and  r^,  since  we  may  write: 

jk  k  ,k  .k 

Q  V  -  -  djrj  - 

In  Appendix  I  we  obtain 


From  the  results  of  Appendices  II  and  III  we  obtain 


N3  ■*  -le^\^(in) 


N4  le^C-Sj - i  K^(lft)) 
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8. 


In  Appendix  III  we  obtain 
a 

^4^“^  ■*  “  "h 

\ 

Furthermore 

/  l^j^(«>)  1^2^"^ 

ft*  I  %T*  M* 

Q  -  Nl  Nj 

Where 

<  -  Jk.c  +  ^(•'k.s  - 

N  and  L  are  similarly  defined. 

a  a 

In  Appendix  I  we  obtain 

•■a  •*•  +  I  ^ 

T*  Sir  /I  1  35  fl  j\ 

Fr  m  ^.!k^  results  Appendices  II  and  III  we  obtain 

N*  >  fle^(K^,(l«)  +  Kj^(iO)) 

2a 
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^2  "*  '  T/2  +  Kj^(if2)) 

2ct 

N*  ^  ne^(Kj,(in)  +  K^(in)) 

In  Appendix  III  we  obtain 

^2^“) 


^^(»)  1 


Let 


Ui  = 


-‘l.S  ■'■  =■'^1.0. 

.3/2 


a 


u<, 


u. 


a 


4.S  i02.c 

a^r2  ■  ■ 

°s  •''  J-°c 
a 


r  ■  ^ — +T!^  (Theodorsen’s  function) 

The  three  equations  become 
“1  •  “2  ■*■  “a  ■  ® 

uj  +  U2  -  -(d*(.  ^  D  +  dj(^)(l.r(l  +  ^)) 
-  ^(2d3  r  -  dj(i-r(i  +  2  ^)) 

3u,  +  U,  -  3u  "0 
12  a 


TECHNICAL  RESEARCH  GROUR 


1 


10. 


Since 


and  r  are  Independent  of  a  and 


-1  1 
1  0 
1  -3 


-8 


u.,,  u  are  all  finite  &8  a  ■*  <»  . 

To  obtain  the  pressure  distribution  on  the  hydrofoil, 

and  also  the  lift  and  moment,  it  is  sufficient  to  obtain  u^  +  U2 

only  since  ocu  5  *»’0  asa-^O  while 
cx  a 


^  . 


;  (see  Appendix  IV) . 


Therefore  the  (normalized)  net  unsteady  pressure  5  on  the  hydro* 
foil  (as  a  ■*  <»)  is  given  by 


^2)  ‘  d*  -  d, 


Similarly  the  (normalized)  unsteady  lift  and  moment  are  given 

by  (see  Appendix  IV) 

/. 


m 


1  ^  *  * 

■  h  "Hc^k  +  ®  V 

•fg 

-►  ir(uj^  +  U2  -  dj) 

-  7r[(^  -  l)d3  -  (^  +  (j^  -  ^)r)dj] 


^1  +  U2  *  dj 

TT  ( - 2 - ®3  ■ 

ir[(^  -l)d3  -(^  +  -  ^)r)dj) 
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To  obtain  the  limit  of  the  fractional  cavity  length 
change  we  write  from  the  first  and  third  of  the  simultaneous 
equations, 


“a - 5 - 

ik  ^ 

Let  0)  ■  +  ItOg  (see  Section  6.) 

ik 

The  0)  -  -  2aa)(a-  +  lo.) 

9  C 

<>  2 

■  -  2a  a)u_ 


2a2 

■  2-”  -  u 
a'^  +  1 


a 


2u. 


Therefore  as  a  •♦  « 


* 

CD 


(Uj^  +  U2) 


w  2  ^ 

For  the  unit  heave  solution  d^  ■  “  10  ,  d^  ■  0 

Therefore: 

f*  -  Tr(2nr  +  in^) 
m*  ■  ir(Or  +  10^) 

(0*  -  -2flr 
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12. 


or 

f*  -  Tr(2nF) 

f*  -  Tr(2nG  + 

m*  ■  irfiF 

m*  -  v(nG  + 

<0*  -  -2fiF 

CO*  -  -2i2G 
c 

where  r  ■  F  +  IG  (F,G  real).  This  Is  Theodorsen' s  result. [2] 

i/t  '  iftf  2 

For  the  unit  pitch  solution,  ■  -  20,  »  -  iO  . 

Therefore 

f*  -  Tr(o  +  or  -  2ir) 

m*  -  Tr(^  +  ^  •  ir  +  ^) 

(0*  -  0  -  or  +  2ir 

or 

f*  -  Tr(n  +  J2F  +  2G) 


<- 

Tr(OG  - 

2F) 

* 

m  " 

8 

G  + 

f) 

* 

m  ■ 

c 

irC^- 

F  + 

f) 

* 

(0  “ 

s 

0  -  OF 

-  2G 

* 

CO  “ 

c 

-  OG  + 

2F 
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This  Is  Theodorsen's  result. [2] 

ilip  4r  ^  ^ 

£  ,  f  ,  m  ,  m  ,  0)-,  0)  were  calculated  In  both  cases 

S  C  8  C  8  C 

(unit  heave  and  unit  pitch)  for  -  0,  .05,  .1,  .3  and  from 
.5  to  4.  in  steps  of  .5  .  The  results  are  Included  in  Table  I 
(under  cu^  -  0). 

Since  to  determine  the  limiting  solution,  the  particu¬ 
lar  term  u^  +  U2  is  sufficient,  rather  than  the  individual 
Uj^,  U2,  u^,  we  see  that  only  the  second  equation  stated  above 
is  required.  This  equation  is  simply  the  statement  that  the 
vertical  component  of  the  velocity  on  the  foil  is  determined 
by  the  foil  motion.  Since  all  physical  models  in  this  context 
Include  this  condition  and  drop  one  of  the  remaining  three 
(pressure  at  infinity,  closure,  continuity  of  the  vertical 
con^onent  of  velocity  in  the  wake),  we  see  that  this  limiting 
behaviour  sheds  no  light  on  the  question  of  the  proper  choice 
of  model. 
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5.  ZERO-FREQUENCT  LIMIT 

We  now  perforio  a  quasl^ateady  analysis  to  determine 
the  behaviour  of  the  solution  in  the  limit  of  low  reduced  frequency. 
This  is  done  by  approximating  the  exact  solution  at  each  time  by 
the  steady  solution  appropriate  to  the  configuration  at  that  time. 
The  solution  obtained  by  this  approximation  is  called  the  quasi- 
steady  solution. 

For  heave  motion  the  quasi- steady  solution  is  essentially 
independent  of  time  so  that  the  zero  frequency  limit  is  zero  for 
all  unsteady  quantities. 

For  pitch  motion^  the  angle  of  attack  varies^  and  a  non¬ 
trivial  result  is  obtained.  If  the  cavity  length  is  less  than  3/4 
of  the  foil  lengthj,  the  magnitude  of  the  unsteady  force;,  moment;, 
and  fractional  cavity  change  are  all  finite.  The  quasi-steady 
nature  of  the  behavior  Implies  that  these  quantities  are  all  in 
phase  (modulo  r)  with  the  motion,  i.e.  the  sine  coii;>onent  is  zero. 
At  3/4  cavity  length  the  magnitudes  of  the  unsteady  terms  are 
infinite,  so  that  a  physically  meaningful  phase  analysis  can  not 
be  made. 

The  steady  state  e]q)ression  for  a  is  given  implicitly  by; 


Therefore 


and 
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The  quasi-steady  change  In  the  length  of  the  cavity 
Is  gotten  from 


^  -  2a>cxAa  ,  where 

i  Ml 

» ui  'i<k+x  M, 

2a  -  /\/?+i 

„  -g  Va^-t-1  !!l 

2a  -  Va^+1 


IT 

o 


Therefore 

60 


JSL 


r2 — 

Va  +1 


h 

M. 


or 


CO- 


(2a  -  Vg  +1 


2a" 


Vg^+i  (2o  -  Vc?+i 


In  the  quasi-steady  limit. 


The  quasi-steady  force  term  Is  derived  as  follows: 


TO  -  ir(7?+r  -  a) 

1  df  _/a  -  Va^+1\ 

TO  ^  7==r-) 

Va^+l 

1  #.  -  Trl^S+LLai 


diU^  ^ 


2a 


-  VTH 
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m 


.(diu2„  )  I-OffSAjl  ^ 
°  2a-V?:i  ° 


or 


*a  -  7r(ya^.fl  +  g) 

ia  -  Vc?+1 


in  the  qiiasl- steady  limit. 


The  quasl'Steady  moment  term  follows  from 


1  <*n  _/  1  -  a 


<aV7 


a  +1  + 


/a  +1  +  a)‘ 
3a  +  2  Va^+1 


4(a^+l)^  (\/a^+l  +  a)^(a^+l) 


777^ 


dm  /..2j.2x_/  3a  +  2\/a^+l  (a^-1)  c/a^+l  +  a)^  \ 

gjj  -  (U  di  )ir(-— a - - 7  + 

(Ja+1  -  2a)(a^+l)  4(Va  +1  -  2a)(a^+l) 


dm  .. 

Am  »  gjj  Ml 


Therefore,  In  the  quasi-steady  limit 

tn  *  -  — == - [a+2(Va^+l+a)  + 

(\/?+r  -  2a)(a^+l)^ 

The  results  associated  with  Q  *  0  In  Table  I 
formulas . 


(a^-l)Va^+l(Va^+l  +  g)^. 
4 

are  based  on  these 
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6.  DISCUSSION 

In  Table  I  we  have  tabulated  the  cosine  and  sine 
components  of  the  normalized  unsteady  force,  moment  and  cavity 

^  *§C  ^ 

length,  f_,  f_,  m^,  m  ,  co_,  oj.,  due  to  unit  heave  and  due  to 

C  0  C  8  C  8 

unit  pitch.  We  may  obtain  the  total  force  f,  moment  m,  and 

cavity  length  cu  from  the  following: 

Force  f  ■  f„  +  Cc(f  cosvt  +  f„  slnvt) 
ore  s 

Moment  m  ■  m^  +  C_(m^  cosvt  +  m„  slnvt) 
o  m'  c  S' 

Cavity  Length  oj  ■  oj^  +  cosvt  +  cOg  slnvt) 

For  heaving  motion  of  amplitude  Bj|^: 

Cf  -  dU^Bj^ 

C„  - 

=0,  - 

For  pitching  motion  of  anq)lltude 

-  diU^Mj 

-  di^U^M, 

m  1 

C  -  (0  M,/M 
CO  o  1'  o 

The  steady  state  textns  are  given  by 
fg  -  -d/U^vCl  + 


-di^U^M„v(- 


a(a  +1) 


4(a^+l)  ^ 


%  ■  l/(a^+l) 
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l<^ 


383? 
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Ok 
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In  these  formulas,  £  is  positive  upwards  perpendicular  to  the 
foil,  m  Is  counter-clockwise  around  the  leading  edge  (with  the 
flow  going  from  left  to  right),  and  co  Is  cavity  length/chord. 

We  now  call  attention  to  some  overall  features  of 
the  data.  In  figures  1,  2,  and  3  are  examples  of  the  type  of 
plot  which  Indicates  the  general  behaviour  of  each  of  the 
unsteady  terms,  constructed  by  plotting  the  sine  coefficient 
("s"  subscript)  against  the  cosine  coefficient  ("c"  subscript) 
and  Indicating  contours  of  constant  frequency  and  constant 
cavity  length. 

All  heave  functions  have  the  following  low  frequency 
behaviour:  amplitude  goes  to  zero,  while  the  phase  tends  to 
t  ir/2,  except  for  the  critical  3/4  cavity  length.  For  this 
case  the  phase  cannot  be  precisely  evaluated,  but  It  appears 
to  be  0  (modulo  v) .  For  hl^  frequencies,  the  force  and  moment 
Increase  In  amplitude  and  tend  to  become  fairly  constant  In 
phase,  the  value  depending  on  the  cavity  length.  The  unsteady 
cavity  length  term  Increases  In  magnitude  at  a  roughly  constant 
phase  up  to  some  critical  frequency  (dependent  on  the  cavity 
length)  and  then  more  slowly  Increases  In  anq>lltude  and  under¬ 
goes  a  rapid  phase  change. 

The  unsteady  pitch  terms  have  significantly  different 
behaviours  at  low  and  high  frequency  (0  *  .5  appears  to  be  the 
transition) .  At  low  frequencies  except  at  the  critical  3/4 
cavity  length,  the  aiiq>lltude  tends  to  the  value  predicted  from 
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the  steady  state  theory,  and  the  phase  tends  to  0  (modulo  tt)  . 
For  the  3/4  cavity  case,  the  amplitude  becomes  oe,  as  expected 
from  steady  theory  while  the  phase  Is  again  hard  to  predict, 
but  seems  to  approach  ±7/2.  At  high  frequencies,  the  behaviour 
of  the  force  and  moment  In  general  resembles  that  of  high  fre¬ 
quency  heave  motion,  l.e.  Increasing  amplitude,  with  no  signifi¬ 
cant  phase  effect. 

The  cavity  length  term  exhibits  an  unusual  behaviour. 
As  we  go  from  zero  to  high  frequency,  for  each  steady  cavity 
length,  the  amplitude  first  decreases  at  roughly  constant  phase, 
then  begins  to  change  phase  while  the  amplitude  reaches  a 
minimum  and  begins  to  Increase.  This  strange  behaviour  may  be 
worth  further  study. 
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7.  CONCLUSION 

Numerical  calculations  have  been  made  to  obtain  the 
quantities  which  describe  the  unsteady  behaviour  of  a  two- 
dimensional  flat-plate  hydrofoil  In  an  Infinite  fluid  without 
gravity  subject  to  oscillatory  heaving  or  pitching,  based  on  a 
particular  physical  model.  The  limiting  behaviour  for  zero 
cavity  length  was  obtained  by  an  asymptotic  development,  while 
a  quasi- steady  analysis  was  carried  out  to  obtain  the  zero- 
frequency  behaviour.  For  zero  cavity  length  we  found  that  the 
limiting  behaviour  does  not  depend  on  the  particular  physical 
model,  l.e.  dropping  the  closure  condition,  that  was  selected. 
The  consistency  of  the  quasi- steady  results  (for  the  zero- 
frequency  limit)  with  the  numerical  results  tends  to  support 
our  choice  of  model. 

It  would  be  of  Interest  to  carry  out  an  asyiiq>totlc 
expansion.  Instead  of  a  quasi-steady  analysis.  In  the  zero- 
frequency  limit  to  determine  If  all  models  are  consistent  In 
this  case  also.  Furthermore  since  experimental  evidence  Indi¬ 
cates  that  the  low  frequency  (reduced  frequency  less  than  .5) 
behaviour  Is  Important,  further  study  of  the  limiting  behaviour 
would  have  Independent  Interest. 

Since  It  has  by  no  means  been  conclusively  demon¬ 
strated  that  this  particular  physical  model  (omitting  the  cavity 
closure  condition)  Is  correct,  although  It  seems  the  most 
reasonable.  It  would  be  of  Interest  to  work  out  all  other  models 
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t^lch  are  possible  In  the  same  context,  l.e.  as  solutions  to 
Laplace's  equation  In  linearized  theory.  Specifically  the  un¬ 
steady  closure  condition  would  be  set  up  as  a  fourth  condition 
and  solutions  would  be  obtained  for  the  models  obtained  by 
dropping  the  continuity  of  the  vertical  component  of  velocity 
and  dropping  the  condition  of  having  the  acceleration  potential 
go  to  zero  at  Infinity.  In  addition  a  third  case  may  be  con¬ 
structed  by  Introducing  an  unsteady  singularity  at  the  trailing 
edge  of  the  foil  (l.e.  dropping  the  unsteady  Kutta  condition) 
and  then  using  all  four  previously  mentioned  boundary  conditions. 
It  would  also  be  desirable  to  extend  the  analysis  to 

Include 

1)  non- zero  thlctoiess  and  camber 

2)  finite  span 

3)  curved,  surface-piercing  foils 

4)  gravity  and  free  surface  effects, 

and  to  perform  a  hydroelastic  analysis  In  each  of  these  cases. 
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APPENDIX  I  -  Asymptotic  Values  of  the  I  Integrals  (see  p.22  and 

^  -Appendix  IV  of  [1]) 

The  Ij^  ^  and  I^  ^  are  of  the  form 

/  cos  <^)  d^ 

-09  r 

2  2  ^  > 
where  P*(4  +a)  +1  (and  where  is  used  instead  of 

for  I*  .) 

’  2 

For  all  values  of  ^  and  a,  P  ^  a  +  1  . 

Therefore  ^  -5S-  ^  so  that  cos(-^)-^l,  for  large  a 

^  a  +1  ^ 

and  fixed  Q,  independent  of  ^  . 

Consequently  the  as3noptotic  value  of  the  integral  may 

be  given  by  the  asyiiq>totic  value  of 

/“  ^  +  .B)  de 

-00  p* 

Similarly  sin  (-^)  ^  ,  so  that 

.2 
P- 


2n  f  a 

^  P 


Since  all  the  Ij^  I|^  g  etc.  are  expressible  in  terms 
of  I^,  its  asymptotic  value  will  first  be  given 
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^  and  Ij^  g  may  be  obtained  from  ^  and  Ij^  ^ 
by  differentiating  with  respect  to  a,  noting  L'Hopltal's  rule. 
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la.c  -  -  3  7 

T  ■»  >  105  tKI 

Ia.s’--53-7 


T  ^  2A  i:idji..g}irjLE).  dg 


2A 


In  3Jn 

(E  -#  +  V) 


Aa'VaN-l 

A 


"•  - ((3a  -  E)  Jn  I, 

rr“  2  1 

((a  +  ^g^+ix^^tl  ■  a).|.  I2a-,:^g^tl} jM 


4a  'Va^+l 


-)Ii 


16a 


L_  (I  +  I 

VTiT  ^  1  ^ 


/  2 

32a(a‘  +  1)  a^  +  1 


8a-^ 
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A.I-5 


Va  +1  ^ 


-  (EI3  +  3J3) 

12  Va +1 

_ AQ _  ,_E  /  3  /  A/?fl+2av 


32  (  ^/a^+l+a) 


•)Ij^+(ori-  ^0^1)33) 


0(ln  e 

•^3 - r  + 1  -^2 


1 

7 


^3,s 


2a-  'Ja^+1  ^  3  ^  Ala^+l+2av 


3a  ^/a^+l  - 


32(/\/a^+l+a) 


+  (a+ ^/c?+I)  -  .^Q^-H+2avv 

32(l+aV 


+1> 


96a  (a+l)‘ 

-  ■C^^a^•f  l-g)  (2a-  Ja^+D+C  \/aKl+a)  (2af 

-1  I  -  r 


1  +  a‘ 


- 7=  (6(-2pl  +  1) 

96aA/a^+l  (o^+lP 


- S3L,  (-2gili_  + _ 1 _ 

32a<a^+l)  ^  ^ 


) 
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A.I-6 


'4,c 


f  de 

*  /"  6^(e^+aKot^-r)  a, 

i-  p5'  o« 


2 

-,3,Q-+2,..  (I  )  (See  Footnote) 
Ha+l) 


-  +  Mj  +  M3 


M, 


2A 


la^/a^+l 


(EI3  +  3J3) 


_  ^  3a*-- 1  ^  _ 1 _ 

64o<a'^  (a  +1) 

A  ^^2v 

-  ( — g—  -  — g-; 


) 


M,  -  — 9 — 

^  2(a>l) 


- -  ((3a^-fF) Jg-FC -^g-t-L-lia) I  ) 

16o(a>l)  ^  4(1+0^)  ^ 


AI1 


V-  (Oa^+FX-fiSj^)  -  F  (,^^^2g)) 
16a(a +1)  ^  4(a +1) 


+  (a-\/oVl  -  2(a*+l))('\/a^+l  -  a  -  li/5L±_2a).)) 

_ *  t  ^ 


a +1 


'  "  ~  ~  ‘T/7  (  -  -»»2  -  ...  . 

128a(?;^ 


Note: 


+  ^a^+1  (3af\Ja^+l  -  2(a^+l)-a^  +  ■?.Cg?tl}r.|g?tjq  ) 

o>l 

(lo  >)  and(l-  _)  stand  for  the  asymptotic  values  of 

OyO 

and  la  -  respectively. 

J,  B 
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128a(a  +1) 


<■ 


jol: 


a  +1  +  a 


+  ( - 2S - 2  + 

A/a  +1  +  a 


a +1 


5 —  3^2  (3a  -  2\la^+l  +  ^  +  3a) 


128a(a  +1) 

M  -  (--Jg  2> 

^  4(a^+l)  8(a'^+l)2 

“l+«2+«3  -  is  f-TT-T  <7^  + 


a) 


”  2^  3/2  “  'n/q^+I)  “  -^Si^SLili 

a(aW^  ./:7:7  ^  (a^+1)^ 


a‘+l 


........1  ■  [3a^-l+(a^-l-l)^-a(a^+l)^/^+3a(a^+l)^/^ 

64a(a  +1) 


(a^+l)^-2a^(3a^+l) ) 


- 2-2 - 5  [^^+2a(a^+l)^/2-2a^(3a^+l)] 

64a(a>l)^ 

- ^ - j  [af(a^+l)^/^-3a^] 

32(a‘^+l)^ 
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I 


A. I 


'4.8" 


N, 


N, 


N, 


dg  +  — A  y  iiu^+sOlsaddi  df 


"  Footnote) 

4(a  +1) 

-  2n(Nj^  +  N2  +  N3) 


2^  EI/+3J 


Va^+l 


8^/7n 


(“IS — > 


(-1—  (T7  I3  -  0^4)  +  3J.) 


.11  E 


8  ^a‘+l  “  a  +1  ^ 


al. 


([- 


8^/c?+l  12  (a +1) 

A  .3aJ3-Fl3^ 
2(a^+l)  ^ 


(ll^;f^)-?3  13.3(3-^)  .3) 


3a^+l 


47^1 


2A  EI3  +  3J3 
( - 12 - ) 


N^l  +  N2  +  N3  -  -  t(E(12o2+ll)-3a(a2+l)-4F(a^+l)3/2 


-4(3a^+l)(a^+l)E)  l3+(9(3(a^+l)3-a(a^+l)E)+12a(a^+l)3/2.i2(3a^+l)  (a^+l))J3] 

"  56(a2+i)5/2  t<E(-12a^-^a^+7)-4F(a^+l)3/2.3o(Qj2^1j)I^ 
+(12a(a^+l)3/2+3(^,2^1)(.3g2^5j.5^(^2^^jgjj^j 

Note:  (13^^) and  (13^5)  stand  for  the  a8ynq>totlc  values  of  I3  ^ 
and  I3  g  respectively. 
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I 

J  A.I-9 

The  coefficient  of 

■  -2a(12aV4a^-7)  +  \|a^+l  (12o^+4o^-7) 

-4a(a^+l)^  +  8^/o^+l  (a^+1)^ 

-3a(a^+l) 

-  -a(28a^-19a^-7)  +  (20a^+20a^+l) 

The  coefficient  of 

-  3(3a^-5)(a^+l)  +  12a  ^^a^+l  (a^+1) 

-18a^(a^+l)  +  9a^|?+l  (a^+1) 

-  -3(a^+lX3a^+5)  +  21a^a^+l  (a^+1) 

N,  +  N,  +  N.  -  - - HTT  [-a(28aVl9a^-7)(2o^+5o^+7) 

^  2  3  I024(a>l)^'^ 

-a(a^+l)  (2a^+4)  (20aS-20a^+l) 

+  '\/c?+l((2aV5a^+7)  (20aV20a^+l)+a^(2a^+4)  (28a^-19o^-7)) 

+  (a^+1) ( 7a(a^+l) ^(6a^+5)+a(a^+l) (6a^+8) (3a^+5) 

-  ^|a^+l  (  (a^+1)  (3a^+5)  (6a^+5)+7a^ (a^+1)  (6a^+8) )  ]  ] 

-  - 1- - -  [Ja^+l  ((a^+l)[20o^(2a^-5a^+7)+2a^(28a^+19a^-7) 

2048  (a  +1)^ 

-  (a^+1)  ( 18a^-45a^+25)  -  7a^  (a^+1)  (6a^+8)  ] 

-  2a^-5a^+7+2a^(28a^+19a^-7))] 

+  a((a^+l)  [7(a^+l)^(6a^+5)+(a^+l)  (6a^+8)  (3a^+5) 

-  (2a^+4)  (20aV20a^+l)  -2a^  (28a^-19a^-  7)  ] 

-  (3a^+7)(28a^-19a^-7)) 
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T 


A.I-10 


[a/oVi  [(a^+l)^[20(2a^-5a^+7)+2(28aS-19a^-7) 


2048(a‘"+l) 


-  ( 18a^+45a^+25)  -  (42aV56a^)  ] 

-  2  (28a^+19a^-  7)  -  (20o^+19)  (2aV5a^+7)  ] 

+  a[(a^+l)^f7(a^+l)  (6a^+5)+(6a^+8)  (3a^+5) 

-  2(20aV20a^+l)-2(28a^+19a^-7)] 

-  2  (a^+1)  (20aV20a^+l)  -  (o^+5)  (28aVl9a^-  7)  j] 

'  2048(0^+1)4  ((a^l)(36o4-3a2-13)-5) 


-  a(  (a^+1)  (-^360^150^+16)  - 136)  ) 


- 2__ -  [(o^+l)( - — - 

2048(0-^+1)^  \  1-2-“ 

(-^  o  +1  -!-  o) 


-  5 


a1c?+1  -I 


I-  136o] 


JL 


2048(0^+1)^ 


[(a^+l)(— -  (o-  \la^+l) 

V?n+o 


3o" 


-  (2o2-(o2+l)-a^/?+I)) 
Ala-^+l  +  o 

-  13a/c?+1  -  16o)  -  5  Vo^+l  +  136o] 

y  i.,,4  ((o2+i)(— Ola^ - ^ 


2048(o"-+l) 


(ori-'\/?+l)2  >/?+!+  o 


+  -  16a)-  5V?;T  +  136  a) 

(^o  +1  +  o)^ 
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2048? 

'"^3 

32a^ 

16? 
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APPENDIX  II  -  Aflymptotlc  VaIum  for  Larea  a  of  J  Integrals  (see 
tl].  P.  17). 

•^k  c  ■  S  /  ^ 

^  •00 

let  w  -  ^ 

o 

then  ^  ■  0  /  co8<^)cD)dw 

*  *00 

similarly 

o 

Jk  a  ■  ^  /  \  8in(nto)dw 

*  •00 

Let 

o 

J  ^  -  n  /  co8(£kD)dw 

Gt «  C  Q 


Jjj  s  -  8in(na>)dw 

^  •  00 


\diere 


^  ^lO  ^  *  °20  -35" 


In  the8e  terms  we  need 


and 


k^l 


for  large  a. 


To  do  this  it  is  first  necessary  to  put  i  in  terms  of  a  and  w 

(see  [1]  page  44) 


4?^  +  4a§^  -  ^  +  I  ■  ® 


2  n - 7 

2€  -  -«+ 


1-1 


a*^l?;rri 

'\j2(a  +  ^a^+1  -  ^ 
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I 


T 


A.II-2 


i>2^  ^2  “  ^  ® 
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I 


A.II-3 


V  1^4  "  ^  (1-  y  (2t^+  ^/?f^  +a)) 

aJo^+I 

7:"^^  9/2  (2(oa/?+I  -30^-1) e2-(2a3+  A/o^+l(2a^+3)) 

4(a  +1)^'  . 


.  ig.qha'J^  -(2a\llah31 
16(a+iy 


TECHNICAL  RESEARCH  9R0UP 


*a)  « 

1) 


4a- 


‘j^O  -  (a^+1)  (of 

uZ  -  I  rT“ 

Qf\j2(cth's}a%l) 


3/2 


^  a 

f20  .  .  3a.^±j±3a.^i!?5l 

^  ^  a^/2(c^f  ^/a^) 

^  .  3al/2‘ 

Therefore 
’^a  a 

(Since,  as  will  be  shown,  a.  and  a  are  0(a),  while  the 

C  o 

1 

e<iuations  involving  ^  and  ^  are  0(1),  •  o(^)  is 

sufficient  to  show  that  these  terms  may  be  neglected.) 

Let  "  *^11  c  ^  '^k  s 


0  / 


iSi» 


dw 


Let  V  ■  -  w 


\  dv 
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T 


Letw^  -  n  /  (/sjl  +  i  -  l)e‘ 


IHv 


dv 


‘«2-^ 


No  •>  w,  a  / 
^  o 


ci^ 


“Ifiv  J 

e  dv 


1  +  Jl  +  - 

A  I  ^ 


a 


00  1  - 


N4  ■*>  W2  a  ^  /  C 


Lllg)  ^  ,-iav 


1  + 


^/^  + V 


dv 


Let  X  -  V  +  1 


Wq  -  fte-^  ^-iftx  ^ 


-  ne^CK^dn)  +  Kj^(in))  +  1 


-  -  ie^^(ifl)  +  (^  +  i) 


'2  ■  ^  (x  -  e"^  dx 

"  T  ^  +  ie^(7  K,,(in)-  ^  K^cin)) 


A.II-5 
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A.  II 


The  results  for  w^,  W2  are  obtained  by  simple  maneuvers 
like  Integration  by  parts,  see  [3]  page  172. 
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! 


A.III-1 


[ 
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A.III-2 


-  - [  (ori-  ^/a^+l)  (2a-  aJc^+I)  -  Oa^+l+SaA/a^+l)  ] 

4a^/o^+l 


- i—  (aVa^P^I  -  1  -  3a^  -  l-3a^/a^+l) 

4a^/cr+l 


1 
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A.III-3 


Part  2: 
B 


•> 


1 

2 


-  B 


2  ^ 


2  /T — 

-y  _  (1  +2+a  \  a  +1 

3  2(a>l) 


1 


'3?  _  2(l-a^)a  A  a^«H  -(2a^-flla^-f3) 

^  16(a^+ir 


•*  - 


^  1  2  39^ 


- 8? -  ((a*+l)  (of  -Jo^+l)  (JSL±2a) 

2aAl?;i . - 


a‘+l 


(3a +l+3a 


-  0 
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A.IV-1 


APPENDIX  IV  -  Asymptotic  Form  of  Pressure.  Lift,  and  Moment  Terms 


A)  Pressure 

Along  the  wetted  surface  of  the  hydrofoil,  the 
pressure  terms  are  functions  of  t)  alone,  where  t)  Is  given  by 


-L 

where  r]  Is  the  value  on  the  upper  surface  and  t]  the  value  on 

-U  • 

the  lower  surface.  Let  A  = 


then 


A 


For  large  values  of  a 


a/o” 

\l^  -  H  > 

^/T" 

1 

a/o” 

Let  6^  ■  ^  »  ^ere  and  are  the  upper  and 

lower  surface  values  respectively. 

Let  6  "  ■  ^  2  C ,  ■  ^  “  ^1.) 

a  ^  ko  da  '^k  ^k' 
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A.IV-2 


1  4.  1 

-  "7+  - 
1  T 


62  "  A 


Observe  that  +  a)  +  1  -  (for  both  signs) 

63  -  +  Wa^+1  -  2o)(ti'^  -  T)‘)) 

A/a +1 


'^3aA 

Therefore 

63  -  .  4  (2^) 


Since  A  ^/a“ 


-  2(-^) 


£ 
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2A 


-  2a)(Ti‘^-Tj")) 


2(at^  +[a^/?fI  -2(a2+l)](T,+-T,-)- 


-  f  (20^)  -  A  (^)(2ah)  -  f  6 


2a 


^lPt-  + 

^  (its  -  (i^)S 


"2 — 7 

xMr-y 

IT 


^  (a^+1)  (art- ^/?+I)  (3a^+l-J-3aVa^+l)  63 


Ba 


<*  VP)-  38“  (- 


Since  B 


2  '\/a 


3  ll~K  ,  3  ll~x 

7a  ^ jh-x  7a  'Vj8+x 


“(i) 
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A.IV-4 


B)  Lift 


Il(-a) 


jA 


I 


1 


R,  «  - 5^ — tj  (2a^+3a^+3+2a(a^+l)^^^) 

3  2(a^+l)^ 

^  2iri 


^ 7  (a(aVl)  (^/a^+l  +  a)  -  4a^) 

^  4(a>l)^ 


(Note:) 


— 

a(a +1) 


77^ 
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A.IV-5 


C)  Moment 


S 


1 


S 


2 


i^(a(a^+13)  H-  JAl  (a^-5))^l 


2(a^+l) 


S 


3 


2 

— M - IT  (4a^+o(a^)+a 'sja^+l  (4aVo(a^)) 

4(a +1)-^ 


X 


64  (a +1) 


(8a^+o(a^)  +  4a  a^+1  (2a^+o(a^))) 


(Note) 


2a(a  +1) 


-gyg  (a(»5a^+19)  +  '\fa^+l  (-5a^+4)) 


Note:  In  (1),  R  and  defined  sligjitly  differently  (i.e.  the 

2 

factor  U  M  was  Included). 
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